Exchange coupling in transition-metal nano-clusters on Cu(OOl) and Cu(lll) surfaces 
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We present results of density-functional calculations on the magnetic properties of Cr, Mn, Fe and 
Co nano-clusters (1 to 9 atoms large) supported on Cu(OOl) and Cu(lll). The inter-atomic exchange 
coupling is found to depend on competing mechanisms, namely ferromagnetic double exchange and 
antiferromagnetic kinetic exchange. Hybridization-induced broadening of the resonances is shown 
to be important for the coupling strength. The cluster shape is found to weaken the coupling via 
a mechanism that comprises the different orientation of the atomic d-orbitals and the strength of 
nearest-neighbour hopping. Especially in Fe clusters, a correlation of binding energy and exchange 
coupling is also revealed. 



I. INTRODUCTION 

The magnetism of transition metal nanostructures on 
metallic surfaces has been studied extensively in the past 
with emphasis on the magnetic moments and ground 
state magnetic configuration. Sophisticated experimen- 
tal techniques for preparation, such as mass-selection 
and soft-landing of free clusters, together with the abil- 
ity to probe the magnetism of these structures on the 
atomic scale, e.g. by X-ray magnetic circular dichroism, 
have considerably advanced the field* 1 "! 2 - Strong fluctu- 
ations of the magnetic properties have been found as 
a function of cluster size and shape, position of the 
individual atoms in the cluster, or geometry of the 
substrate ^^^' 7 ' 8 ' 9 ' 10 ' 11 ' 12 In some cases, though, it 
has been possible to recognize and interpret regularities, 
and derive rules of thumb based on an understanding of 
the electronic structure. 

Supported clusters of atoms from the end of the 3d- 
series (Fe, Co, Ni) are known to be ferromagnetic, due to 
a double-exchange mechanism. The spin moments in Fe 
clusters have been found to vary in a regular way, namely 
the moment of each Fe atom scales down linearly as a 
function of the number of Fe neighbours! 8 ' 9 ' 10 This nice 
effect seems to be independent of the substrate, and per- 
sists also in FeCo clusters^ On the other hand, in pure 
Co clusters the moment is practically saturated, while in 
Ni clusters the atomic moments do not seem to correlate 
with the coordination number. 

For V, Cr and Mn clusters, it has been found that 
the intra-cluster coupling is antiferromagnetic, leading 
to frustration and non-collinear magnetic order if the 
substrate provides appropriate geometry (e.g., (Ill) 
surfaces) ! 13 ' 14 However, if the substrate is magnetic, the 
exchange coupling to the substrate atoms also plays a role 
and non-collincar magnetism can appear in Cr and Mn 
clusters or chains on (001) Ni and Fe surfaces ! 15 ' 16 ' 17 ' 18 

The magnitude of the exchange coupling in 3d clusters 
on surfaces has been less studied so-far. This is an im- 
portant quantity, however, in determining the stability 
of the ground state and the crossover or blocking tem- 
perature. In addition, recent developments in experi- 
mental techniques allow to probe the spectrum of mag- 



netic excitations, e.g., via inelastic scanning tunneling 



spectroscopy 
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or spin-polarized scanning tunnel- 



ing spectroscopy* 

From the theory point of view, calculations on 
supporte d 10 ' 24 and free-standing 2 ^ transition-metal clus- 
ters have shown that the inter-atomic exchange coupling 
fluctuates strongly with respect of cluster size, shape, or 
position of the atoms in the cluster. These works moti- 
vate an analysis of the local electronic structure effects, 
in order to obtain a better insight into the driving mech- 
anisms for these fluctuations. 

In this paper we perform such an analysis. We present 
density-functional results on Cr, Mn, Fe, and Co nano- 
clusters (up to 9 atoms large) of mono-atomic height on 
Cu surfaces. For the ferromagnetic clusters (Fe and Co) 
we also consider Cu(lll) as substrate, while for the an- 
tiferromagnetic ones we restrict our study to collinear 
magnetism on Cu(001) (Cr and Mn show non-collincar 
magnetism on Cu(lll)). Our focus is on the inter-atomic 
exchange coupling, where we attempt to find and under- 
stand the trends via the details of the electronic struc- 
ture. It is gratifying that, in some simple geometries, the 
results can be interpreted in a transparent way within a 
simple tight-binding model. 



II. METHOD OF CALCULATION 

Our calculations arc based on density-functional the- 
ory within the local density approximation (LDA), as 
parametrized by Vosko, Wilk, and Nussair^i The Kohn- 
Sham equations are solved in the framework of the 
full-potential Korringa-Kohn-Rostoker Green-function 
method (KKR) with exact treatment of the atomic cell 
shapes^ 7 - using an angular momentum cutoff of Z max = 3. 
The Cu surfaces were modelled by slabs of a finite thick- 
ness of 18 atomic layers; in all calculations, the LDA equi- 
librium Cu lattice parameter (3.51 A) was used, while 
structural relaxations were neglected. Within the KKR 
method, first the Green function of the host (Cu surface) 
is calculated, while in a second step a Dyson equation 
is solved for the Green function of the embedded cluster. 
For all clusters, neighbouring host atoms were considered 
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in the self-consistent calculation in order to account for 
the screening of the charge. 

The exchange coupling is calculated by the method 
of infinitesimal rotations^ This method is based on a 
hypothesis of correspondence between the energy change 
AH of the Heisenberg Hamiltonian 
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and the total energy change found within DFT, AEdftj 
upon rotating the magnetic moment directions e n and 
e n '. By virtue of the magnetic force theorem this 
energy change can be calculated by the difference in 
the Kohn-Sham eigenvalues. One obtains the ex- 
change constants by taking a second derivative, J nn i = 
-\d 2 H/(de n de n ,) = -±d 2 E BFT /{de n de n ,), as^ 
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where Gli (E) is the inter-site KKR structural Green 
function for spin up (f) or down (J.), are spin- 
dependent scattering matrices at sites n and n', and Tr; m 
indicates a trace in angular momentum indices (G nn > 
and t n are matrices in angular- momentum space). We 
have also identified the integrand as an "exchange cou- 
pling density" j n n'(E), which is useful for the analysis 
of the results. In short, if the moments of two atoms 
arc rotated with respect to each other by a small angle 
9, then the difference in energy arises primarily due to 
shifts of the atomic energy levels or resonances. Then, 
2 j nn '(E) gives the energy-shift of the local states at E, 
while 9 2 J nn i gives the change in total energy; a subtrac- 
tion of the interaction with the rest of the magnetic atoms 
is also included implicitly. We also introduce the sum of 
the exchange interactions in the n-th atom in the cluster, 



r(«) 



E 



Jnn' i 



(4) 



which is a measure of "local spin stiffness" . Integrating 
j up to a certain energy E, instead of Ep, yields J(E) 
and Jo (E) that also helps to analyze the contribution of 
specific states. 

In practice, the set of equations ([1]) and ([2]) is an 
accurate parametrization of the DFT total energy only 
for small deviations from a reference state. Usually the 
ground state is chosen as reference, which can be ferro- 
magnetic or antiferromagnctic. Eq. ([2]) then yields J > 
if energy must be payed for mutually rotating the spins 
direction, and J < if energy is gained. Thus, starting 
from a ferromagnetic or antiferromagnetic state within 
the LDA, J > always indicates stability of this state, 
irrespectively if it is ferro- or antiferromagnetic; the ex- 
tra (— ) sign associated with antifcrromagnctism in JTJ) is 
absorbed, so to say, in either e„ or e„/ . 
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FIG. 1: Atomic spin moment of Cr, Mn, Fe, and Co in corre- 
sponding clusters of 14 different cluster geometry types, rang- 
ing from 1 to 9 atoms, on the Cu(001) surface. The moment 
shows a linear trend as a function of the coordination number 
of nearest transition-metal (TM) neighbours. 



III. ATOMIC SPIN MOMENTS 

In Fig. [1] we show the dependence of the atomic mo- 
ments as a function of neighbouring atoms for Cr, Mn, 
and Fe, and Co. The strongest atomic spin moments are 
found, as expected, in Mn clusters, ranging from 4.23 fiB 
for a single adatom to 2.89 /ie for the central atom of 
a 9-atom cluster. Cr shows somewhat lower moments, 
then comes Fe and finally Co. In a previous paper— we 
showed that the atomic moments M in Fe drop with in- 
creasing coordination following approximately a simple 
linear relation, 



M(i) = -aN c (i) + b, 



(5) 



where N c is the number of nearest neighbours of the i- 
th Fe atom in the cluster and a and b are positive con- 
stants (independent of i). This relation simply expresses 
the common wisdom that the moment decreases with hy- 
bridization, and that the hybridization mainly depends 
on the number of nearest neighbours. However, this was 
found not to be so simple for Co, where the moment is 
more saturated and for Ni4 Additional effects come from 
interference among the d states of remote atoms in the 
cluster, at the onset of a narrow d-band formation at the 
Fermi level Ep, giving a stronger relative scatter to the 
data. 

In the case of Cr and Mn clusters, relation |5]| holds 
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rather well (here for the absolute value of M, since the 
clusters are antiferromagnetic). As the d-states of these 
elements are either below the Fermi level (for majority- 
spin) or above the Fermi level (for minority-spin), there 
is no e?-band formation at Ep (see Fig.[2]for the density of 
states of Cr, Mn, Fe, and Co dimers). In addition, the an- 
tiferromagnetic configuration allows for d-d hybridization 
only between occupied and unoccupied states making it 
comparatively weak, which is reflected on the smooth, 
Lorenzian-like, minority-spin density of states, contrary 
to Fe and Co. (The majority-spin resonances appear 
wider in Mn, Fe, and Co, because of their interaction 
with the Cu substrate d-band that lies at approximately 
the same energy.) Since the hybridization at Ep arises 
primarily due to the extended s states of the first neigh- 
bours, the behavior of the moments correlates smoothly 
with the coordination. In the case of Co, the moment 
is almost saturated, changing very little with coordina- 
tion; the fluctuation in M cannot be correlated to N c . 
The weak scatter of data for any particular coordination 
corresponds to different atoms on different clusters. 



IV. EXCHANGE COUPLING: 
REMARKS 



PRELIMINARY 



We proceed with the exchange coupling in the nano- 
clusters. Before presenting the density-functional results, 
we discuss how one can interpret these in terms of a sim- 
ple model. 

The inter-atomic exchange coupling depends on vari- 
ous factors. The nearest-neighbour exchange can be in- 
terpreted in terms of the local density of states (DOS). 
The longer-range interactions, on the other hand, are 
governed by a generalization of the Rudcrman-Kittcl- 
Kasuya-Yosida (RKKY) interaction,— depend on the 
properties of the Fermi surface, and are far more complex 
in nano-clusters, as a Fermi surface is not yet formed. 
Also, details of the surroundings of the cluster can play 
a role due to quantum confinement.— Furthermore, the 
substrate below the surface layer, e.g. buried magnetic 
clusters or layers^ can affect the long-range interactions; 
even non-magnetic buried structures can play a role as 
spin-dependent wavefunctions incident from a magnetic 
atom at the surface can be reflected at the buried struc- 
ture and propagate to another atom back on the surface, 
especially in the presence of Fermi-surface focusing^ 
Another type of exchange interaction is the anisotropic 
exchange of the Dzyaloshinskii-Moriya type, which has 
been found, for example, to produce a long- wavelength, 
non-collincar magnetic ground state in Mn overlayers on 
W(110)i^ This, however, depends on the spin-orbit cou- 
pling, and even in transition-metal clusters on Au and Pt 
(where the spin-orbit coupling is strong), the anisotropic 
exchange is found to be two orders of magnitude weaker 
than the nearest-neighbour Hcisenberg exchange. 34 ' 35 In 
Cu, spin orbit coupling is much weaker, therefore we con- 
sider that its effect can be safely neglected for the small- 
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FIG. 2: Density of states and exchange coupling for Cr, Mn, 
Fe, and Co dimers on Cu (001). Label "Fe-de" corresponds 
to a Fe dimer with the majority states pushed to low en- 
ergies, so that the double exchange becomes evident. Full, 
black line over shaded area (left scale): density of states; 
upper panels correspond to majority spin, lower panels to 
minority spin. Dashed, blue line (right scale): exchange cou- 
pling J(E) — J j(E')dE' as a function of energy; J(Ef) 
corresponds to the actual exchange constant. Arrows indi- 
cate the particular local maxima of J(E) that correspond 
to the kinetic-exchange mechanism (absent in the "Fe-de" 
case, where both maxima of J(E) correspond to the double- 
exchange mechanism) . Dashed-dotted red line: exchange den- 
sity j(E) (not to scale). 
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size clusters that we examine here. Following these con- 
siderations, we focus our discussion only on the nearest- 
neighbour Hciscnberg-type exchange coupling, which is 
also much stronger than more-distant interactions. 

As is long known^ there are two main contributions to 
J, of different origin: the first, sometimes called "kinetic 
exchange" , favours antifcrromagnctism in Cr and Mn. 
The second, sometimes called "double exchange" , favours 
ferromagnetism in Co and Fc. Both mechanisms are asso- 
ciated with direct hopping among d states of neighbour- 
ing atoms. As we will see below, the two mechanisms can 
be present simultaneously and compete, reducing the ab- 
solute value of J.— 

To proceed with the discussion, consider two neigh- 
bouring transition metal atoms. Kinetic exchange arises 
from a level repulsion between occupied majority states 
of an atom with unoccupied minority states of a neigh- 
bouring atom, when the moments of the two are op- 
positely oriented. In a simple tight-binding picture, let 
us assume that the energy of the occupied d level is at 
Ed — A < 0, while the unoccupied level is at Ed + A > 0, 
where 2A is the exchange splitting; the Fermi energy 
is assumed to be at Ep = 0. Allowing for a hop- 
ping t between two neighbouring atoms, the levels move 
to E± = Ed ± yj A 2 + t 2 . Upon forcing the two mo- 
ments to a ferromagnetic configuration, the shift van- 
ishes and the occupied levels move higher, which costs 
energy 2[(E d - A) - E-] = 2(VA 2 + t 2 - A) « t 2 /A, for 
A 3> t (the factor 2 accounts for the number of occupied 
states). This picture does not change even in the pres- 
ence of a broadening by a hybridization T of the levels 
with a background continuum (as long as the broadened 
resonance does not cross Ep), since the center-of-mass of 
the d levels is normally not shifted by T. 

Double exchange arises if the majority-spin or 
minority-spin states of the two transition metal atoms 
are in the proximity of Ep, i.e., either Ed — A ks or 
Ed + A w 0. Suppose that the two moments are ferro- 
magnetically oriented, and that Ed + A = 0. Allowing 
for a hopping t between the d-states, the minority lev- 
els split into bonding and antibonding, and are driven to 
E± = ±t; E- is occupied, E+ empty. Upon forcing the 
two moments to an antifcrromagnctic configuration, the 
splitting vanishes, and the occupied EL move back to Ep. 
This costs energy t. (The same splitting occurs also for 
the majority levels, but does not affect the total energy 
as they are fully occupied and their center-of-mass does 
not change). The energy gained by the double-exchange 
mechanism, however, is reduced if the d-states are res- 
onant with a background hybridization T, because the 
tails of the bonding and antibonding resonances cross 
Ep, so that their rcpopulation partly counter-acts the 
energy gain. A simple calculation shows that the energy 
gain by the hopping, J^ oo EAn(E)dE, is a decreasing 
function of T (An(E) is here the difference between the 
density of states without and with hopping, and only the 
single-particle energies are considered). 

Now, if the double-exchange mechanism is present, it 



must coexist with kinetic exchange, which cannot be 
switched off. There is then a competition between the 
two, and in the simple picture given above (disregard- 
ing the hybridization), the criterion for a ferromagnetic 
ground state is t > t 2 /A, i.e. A > t. Otherwise the 
ground state is antiferromagnctic. This analysis also sug- 
gests a way to reveal the strength of the double-exchange 
mechanism in a calculation, which we use below. If the 
majority-spin states are artificially driven to lower ener- 
gies by acting on them with an attractive potential, the 
kinetic exchange t 2 / A will decrease and what is left will 
be mainly the double exchange part. 

To conclude this subsection, nearest-neighbour ex- 
change coupling can partly be understood in terms of 
bond formation, but there are competing factors that 
decide its final value. 



V. EXCHANGE COUPLING: RESULTS AND 
DISCUSSION 

We now present and analyze the exchange coupling 
calculated within density-functional theory. Results on 
most of the investigated clusters are shown in Fig. [3] 
As we expected, we find an antiferromagnetic coupling 
for Cr and Mn clusters, dominated by kinetic exchange, 
while Fe and Co are ferromagnetic, dominated by double 
exchange. On the average, the interaction strength de- 
creases with increasing cluster size, dimers showing the 
strongest coupling, but there are many exceptions (es- 
pecially for Mn), while the local geometry around each 
atom plays a role. 

Cr is the champion in the strength of exchange cou- 
pling, with values as high as |J| = 134 meV for the 
Cr/Cu(001) dimer. Mn clusters show comparatively 
weak coupling: for the Mn/Cu(001) dimer we find | J| = 
28 meV, a value which increases by up to 70% for certain 
configurations. In Fe the coupling is again strong, reach- 
ing J = 81 meV for the Fe/Cu(001) and 100 meV for the 
Fe/Cu(lll) dimer. In Co, finally, the coupling is again 
relatively weaker, reaching values of J = 68 meV for 
the Co/Cu(001) and 80 meV for the Co/Cu(lll) dimer. 
We proceed with an analysis of our results based on the 
model that were presented in the preceding section. 



A. Dimers: exchange coupling and density of states 

In a first step, we discuss the exchange coupling in 
dimers on Cu(001) and relate it to the electronic struc- 
ture. The picture that is derived here is helpful in un- 
derstanding the behavior of larger clusters. 

Our guide for the discussion is Fig. O There, the the 
density of states is shown for one of two atoms of a Cr, 
Mn, Fe, and Co dimer (the other atom has exactly the 
same DOS, but for Cr and Mn the spin directions are 
interchanged). In the same figure, the energy-dependent 
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Fe clusters on Cu(111) 
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FIG. 3: Nearest-neighbour exchange coupling (absolute value, in meV) for Cr, Mn, Co, and Fe nano-clusters on Cu(001) and 
Cu(lll). A schematic top- view of the clusters is shown, the Cu surface is understood to be one atomic layer below. The view 
corresponds to surface- adapted geometry, i.e., rotated with respect to the fee cubic axes by 45°. Coupling values identical by 
symmetry are not annotated more than once, but indicated by the same line-type (e.g., full or dashed line). For Co/Cu(001), 
Cr/Cu(001), and Fe and Co/Cu(lll), only selected clusters are shown. On the (111) surface, couplings that appear on first 
sight equivalent by symmetry can differ, due to the relative displacement of the sub-surface layer. 
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exchange coupling density. j(E), and coupling, J(E) = 
J E j(E')dE', are shown (see Eq.[3J). 

In all cases, the majority- and minority-spin d- 
resonances can be clearly seen. Around each resonance, 
j(E) (dash-dotted red line) shows an S-shapcd form, re- 
vealing the change of the resonance form upon rotating 
the spins with respect to each other. In Cr and Mn, 
the lower levels of each resonance will move even lower 
and the higher levels will move higher: the resonance will 
widen up. This is because the starting point is here an 
antifcrromagnetic orientation of the moments, and rota- 
tion allows for majority-majority and minority-minority 
hopping. In Fe and Co, where the starting point is a 
ferromagnetic orientation, the opposite will happen: the 
resonance will narrow down, because majority-majority 
and minority-minority begins to be blocked; thus, the 
S-shape is inverted. As the minority spin resonance is 
bisected by E F for Fe and Co, the single-particle part of 
the energy, / F E n{E) dE, will increase, signalling a fer- 
romagnetic coupling of the double-exchange type. Note 
that the behaviour of j(E) is more spiked at the majority 
resonance in Mn and Fe, due to the hybridization with 
the Cu d-states; moreover, in Cr, Mn, and Fe the shape 
of j(E) around the minority-spin resonance appears to 
have a double- S form, revealing contributions from dif- 
ferent d-orbitals at slightly different energies. Around the 
S-shapes, the integral J(E) shows peaks (positive or neg- 
ative, depending on the starting moment configuration) 
revealing the total strength of the interaction. 

At some energy between the resonances, j(E) passes 
through zero; however, J(E) can show a plateau with a 
local maximum at those energies. This reveals a net shift 
of the resonances, apart from the broadening or narrow- 
ing, upon rotation of the moments. The net shift stems 
from the majority-minority-state repulsion, and can be 
identified with the kinetic exchange. Its maximal value, 
indicated by arrows in Fig. [21 is strong in Cr, but less so in 
Mn, probably because the Mn d-states are lower in energy 
and thus more localized spatially, reducing the hopping. 
This is one reason why Mn shows a weaker exchange cou- 
pling than Cr. Another, more important reason is that 
for Cr E-p bisects J(E) at a high point, in the middle of 
the plateau, while for Mn E-p is already at a point strong 
descent of J(E), close to changing sign. Here, the double- 
exchange mechanism among the minority-spin states is 
already starting to set in, competing with the kinetic 
exchange.— This affects the behaviour of the coupling 
with increasing cluster size, as we discuss later. 

The kinetic exchange is also present for Fe and Co; it 
is just that double-exchange is stronger and dominates, 
causing a ferromagnetic ground state. It is interesting to 
disclose the strength of the double-exchange by shifting 
the majority-spin resonances of Fe and Co to lower en- 
ergies, so that the kinetic exchange is diminished. We 
achieved this by acting on these states with an attractive 
potential of 2.7 eV. This resulted in an increase of the 
coupling by approximately 50% for all Fe and Co clus- 
ters. An energy-resolved picture of this can be seen in the 



case of a Fe dimer, the DOS of which is shown in Fig. [2j 
labelled "Fc-dc" . Here, J(E) docs not show a peak be- 
tween the resonances any more (the kinetic exchange is 
now absent), while it is bisected by E-p much closer to its 
maximum than in the normal Fe dimer. 



B. Trimers: the role of local geometry 

Having discussed the exchange in dimers, we now focus 
on the complications brought about by local geometry 
effects. These can be best demonstrated in the example 
of trimers, since for larger clusters there are too many 
parameters that play a role. 

We begin with an observation. The exchange coupling 
in linear trimers is weaker than in corner-shaped ones 
(see Fig. EJ). This is true for Cr/Cu(001), Fe and Co on 
Cu(001) and Cu(lll), but not for Mn/Cu(001). This 
behavior can be interpreted again in terms of a tight- 
binding model, but with more than one orbitals. We will 
do so, deferring the explanation of the Mn to a later point 
in this subsection. 

In the interaction causing the kinetic or double ex- 
change, not all d orbitals take part equally. Some can 
have an orientation not favouring bonding, while oth- 
ers can be pointing directly into the neighbouring atom. 
Consider for instance a linear-shaped trimcr in the y- 
direction, setting for definiteness the x and y axes in 
the surface plane (schematic positioning of atoms 1, 2, 
3 in Fig. Hfc). The d xy , d yz , and d x 2_ y 2 orbitals point 
partly toward the neighbour atoms and can hybridize 
with their counterparts there. The d z 2 orbitals point out- 
of-plane, and can therefore hybridize only weakly with 
d x 2_ y 2 neighbouring orbitals (hybridization with the rest 
is forbidden by symmetry, if one disregards the surface- 
induced symmetry breaking in a first approximation). Fi- 
nally, the d xz orbitals point away from the neighbours, 
and are practically non-bonding. Thus three channels 
have a major contribution to the exchange coupling: one 
consists of the three d xy orbitals, one stems from the 
three d yz , and one from the three d x 2_ y 2. Because each 
of these channels is brought about by three orbitals, we 
call them "triple channels" henceforth. If the trimer were 
to be extended to a linear chain, these three triple chan- 
nels would form the most dispersive d bands. 

Now let us change the picture to a corner-shaped 
trimer (schematic positioning of atoms 1, 2, 3 in Fig. SJi). 
Here we have only two triple channels, namely the d x 2_ y 2 
and the d xy channels where all atoms participate. How- 
ever, we also have two bonds: d yz -d yz between atoms 
1 and 2, and d xz -d xz between atoms 2 and 3. In effect, 
compared to the linear trimer, one triple channel has been 
replaced by two bonds. 

A straightforward calculation shows that the energy 
gain of two simple bonds is larger than the one of a triple 
channel. Considering first the case of double exchange, 
assume that the majority-spin d levels of the three atoms 
lie at Ed + A = Ep = 0. In the linear trimer, if we open 
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a triple channel through the d yz states by allowing for a 
hopping t, a threefold splitting is created with E± = \/2t, 
Eq = 0; the d xz orbitals are not affected. At half-filling 
of the minority-spin orbitals (optimal case for double- 
exchange), the energy gain is AE^ ceiT — \/2t. In the 
corner-shaped trimer, allowing a hopping between d yz - 
d yz of atoms 1 and 2 and d xz -d xz of atoms 2 and 3 creates 
two double splittings of E± = t. Now, at half-filling of the 



minority-spin orbitals, the energy gain is A£^° rner = It. 
Evidently, A#£° rner > A££j t nGar . Carrying out the cal- 
culation for the kinetic- exchange case at full majority- 
spin and empty minority-spin orbitals gives AE^ ea,T = 
2(VA 2 + i 2 - A), AEgf™ = 4(VA2 + t 2 - A), i.e., 
again AEg° rnoY > A££j t near . In both cases, the linear 
trimer should show a smaller exchange J than the corner- 
shaped one. Note that no second-neighbour hopping is 
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necessary in this interpretation. On the (111) surface, 
the corner-shaped trimer is not in a cj> = 90° configura- 
tion, but rather at <p = 60° or <j> = 120°. Then the effect 
is still present, but to a lesser extent as the hoppings sat- 
isfy partially a triple channel and partially a bond (the 
hopping contributing to the triple channel is t cos <fi, the 
one contributing to the bond is tsmcj)). 

This simple picture can be recognized in the m-resolved 
density of o?-states. In Fig.|3]we show the minority-spin d- 
DOS for a Fe/Cu(001) adatom, dimer, linear and corner- 
shaped trimer. In panel (a), among the resonances of 
the single adatom, d xy and d x 2_ y 2 are clearly the most 
localized, as the lobes of these orbitals extend mainly in- 
plane. Upon adding a neighbouring Fe atom (panel b), 
clear-cut splittings arc formed for all orbitals that point 
toward the Fe neighbour. Even the d z ? orbital hybridizes 
somewhat with the d x 2_ y 2 of the neighbour, causing a 
three-peak structure in d x 2_ y 2. Only the d xz orbital re- 
mains oblivious to the presence of the second atom, as 
it points in the wrong direction and is orthogonal to any 
neighbouring orbital pointing in its way. 

Next we add a third atom to form a linear trimer (up- 
per and lower panel c). Here, more splittings appear. In 
the two outmost atoms, 1 and 3, orbital d xy shows now a 
threefold splitting, corresponding to the three levels E± 
and Eq that are formed within a triple channel, as we dis- 
cussed above; the same is true for d yz . This is the onset of 
band formation. However, the state corresponding to Eq 
has no weight in the middle atom, which can also be ver- 
ified by solving the associated 3x3 system of equations 
within the tight-binding model. For the d x 2_ y 2 orbitals 
the picture is more complicated due to their additional 
hybridization with d z 2 . 

Finally, we move one atom to form a corner-shaped 
trimer (panel d). Here the d yz triple channel is broken. 
Instead, we have a bond between the d xy of atoms 1 and 

2, and an identical bond between the d xz of atoms 2 and 

3. The splittings, looking exactly like the d yz splitting 
of the dimer of panel (b), verify this picture. We con- 
clude that the tight-binding model indeed interprets the 
calculated DOS and coupling in Fe and Co. 

We turn now to the interpretation of the behavior of 
Cr and Mn. For Cr, which is a kinetic-exchange system, 
it is not the level splittings but the level repulsion that 
causes the 90°, corner-shaped configuration to show a 
stronger coupling than the linear configuration, as out- 
lined above. Mn, on the other hand, where the compe- 
tition of kinetic and double exchange is important, be- 
haves differently. Upon going from the linear trimer to 
the corner-shaped one, both the kinetic and the double 
exchange are strengthened. Apparently, however, the en- 
hancement of the latter is stronger, so that, in the end, 
they sum up to a weaker antiferromagnetic coupling. 

To conclude this subsection, a simple, nearest- 
neighbour tight-binding model can be used to qualita- 
tively interpret many of the features of the DOS in dimcrs 
and trimers. This supports the interpretation of the 
exchange coupling that we developed earlier, without 



a need to introduce second- neighbour corrections; such 
corrections are of course always present, but apparently 
not of central importance. The complications of local 
geometry effects are revealed and explained in the case 
of trimers; in particular it becomes clear that very im- 
portant for the behavior is not only the homogeneous 
broadening of the resonances, but also the onset of band 
formation creating geometry-dependent splittings. If the 
moments are tilted away from parallel orientation, these 
energy-gaining splittings will start closing up, giving rise 
to the exchange coupling. We now proceed to the discus- 
sion of the exchange in larger clusters. 



C. Increasing size: effect of hybridization 

Except for the case of Mn, increasing cluster size 
leads to weaker coupling. This effect is coming from 
the stronger hybridization with increasing coordination. 
For Fe and Co, hybridization contributes to the lowering 
of exchange coupling in two ways. First, the double- 
exchange interaction is weakened, as we described in sec- 
tion IIVI Second, the tendency for magnetism is weak- 
ened, as reflected also in the reduction of the atomic mo- 
ments of Fe. 

For Cr clusters, which arc a kinetic- exchange systems, 
mainly the latter effect should play a role, since kinetic 
exchange should depend on the splitting but not on the 
hybridization, as we argued in section US However, for 
this argument it was hypothesized that the resonance 
does not cross E-p, which is not completely true here. 
As one can see from the density of states (Fig. the 
resonances do cross E-p to an extent that increases with 
hybridization (reflected also in the decrease of the Cr 
moments). Thus the repopulation of the resonance tails 
also contributes to the weakening of the kinetic exchange. 

The effect of hybridization becomes clearer if one com- 
pares the values of the exchange coupling on Cu(001) to 
the ones on Cu(lll). In the former case, each atom in the 
cluster has four Cu neighbours, in the latter only three. 
We see that, on the average, the coupling on Cu(001) 
is indeed weaker than on Cu(lll) (for the same number 
of transition-metal neighbours). Further evidence comes 
from examining Fe and Co clusters on Ni(001) and (111) 
surfaces (not shown here in detail). Here, additional 
background hybridization is induced on the one hand 
from the Ni d states at Ep and on the other hand from the 
smaller lattice parameter. Although the moments, as- 
sisted by the ferromagnetic substrate, are slightly larger 
than on Cuf£ the exchange coupling within the cluster 
is significantly weaker, by a factor of up to 50% for 
Fe/Ni(001) compared to Fe/Cu(001), 40% for Fe/Ni(lll) 
compared to Fe/Cu(lll), and 25% for Co/Ni(001) com- 
pared to Co/Cu(001). 

It is also interesting to compare the systems 
Co/Pt(lll) to Co/Au(lll), calculated by Sipr et alJ£ 
In order to eliminate the effect of increased lattice pa- 
rameter, Sipr at al. also calculated Co/Au(lll) in the 
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lattice parameter of Pt. Although the Pt substrate is po- 
larizable, enhancing the moments of Co/Pt compared to 
Co/Au, the pair exchange coupling is weaker for Co/Pt 
dimers and trimers, compared to the same formations on 
Aui^ This effect could be due to the strong hybridization 
of the Co d states with the Pt d bands at Ef (the d bands 
of Au are considerably lower). For larger clusters, how- 
ever, this trend is not so clear and depends on the local 
geometry. Additional evidence on the effect of hybridiza- 
tion comes from calculations of Co ad-clusters on Pt(lll) 
in Ref. [24] where it is found that the nearest-neighbour 
exchange of Co atoms at the cluster surface are larger by 
approximately a factor 2 compared to interactions of a 
central Co atom. 

Finally we comment on the irregular behaviour of Mn 
clusters. The irregularity arises from the competition 
of the kinetic and double exchange analyzed in the dis- 
cussion on the dimers above, and we also saw it in the 
case of trimers. Small changes in the electronic struc- 
ture, caused e.g. by local geometry effects or cluster size, 
can shift the turning point of J(E) up and down in en- 
ergy, disturbing the delicate balance between kinetic and 
double exchange. This is the reason that the Mn clus- 
ters show a rather irregular behavior in J nn >. More to 
this point, increasing coordination and hybridization will 
suppress the double-exchange mechanism stronger than 
the kinetic exchange, as we discussed in Sec. IIVI There- 
fore, the exchange coupling tends to increase with size in 
Mn clusters, opposite to Cr, Fe, and Co. 



D. Correlation of exchange and formation energy 

As we discussed so-far, level splittings and shifts are 
responsible for the nearest-neighbour exchange coupling. 
However, it is also known that such splittings are con- 
tributing to the binding energy. Particularly for d-states, 
this idea dates back at least to the model of Friedcl for 
the cohesive energy of transition metals. Therefore, there 
could be a correlation between the energy gain upon at- 
taching an atom to a cluster and the exchange coupling 
of this atom with its neighbours. However, other factors 
also enter the formation energy, and they could overweigh 
the effect of the d states. 

Of course there is a trivial correlation, in the sense that 
the total exchange coupling energy of an atom increases 
with the number of neighbours, and so does its binding 
energy. What we mean here is a correlation even for 
atoms of the same coordination. To be more precise, we 
are looking for a linear relation, of the form 



E b {n) 



E r , 



(6) 



where E b (n) is the binding energy of the n-th atom to 
the rest of the cluster and A and i? res t are constants. 

In fact we do find such a correlation for Fe/Cu(001) 
and Fe/Cu(lll), but not for Cr, Mn, or Co (apart from 
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the trivial one mentioned above). We calculate the bind- 
ing energy E b = E N+X - (E N + E ad ), where E N is the 
total energy of an TV-atom cluster, E^+i the total energy 
of the same cluster augmented by one atom, and £" a d the 
total energy of a lone-standing adatom (total energies 
are calculated with respect to the reference energy of the 
clean Cu surface). 

In Fig. [5] we show the relation between the absolute 
value of the binding energy of cluster atoms, |-Eb(n)|, and 
the total exchange coupling of these atoms to the rest of 
the cluster, Jfi n ' = ^ , J nn i . We only show results for 
atoms at the cluster edge, as we consider it unlikely that 
a cluster will be formed with an atom missing in the 
middle; therefore the maximal number of neighbours is 
3. We see a clustering of the data according to the num- 
ber of neighbours of the added atom, as expected, but 
within each of these groups, a correlation of the form © 
is evident. There is also a clear "background" coupling, 
termed -E res t in relation which varies with the num- 
ber of neighbours. Deviations from relation © can be 
caused by several effects, such as participation of the s 
electrons in the bonding, charge relaxation, moment re- 
laxation (which is rather strong for Cr and Mn), etc. In 
addition, in the case of Mn, Fe and Co, double exchange 
and kinetic exchange compete in the final value of J nn i , 
while the level shifts associated with kinetic exchange do 
not contribute to the binding energy for Fe and Co, and 
the shifts associated with double exchange do not con- 
tribute to binding in Mn. 
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VI. CONCLUSIONS 

The exchange coupling in supported transition metal 
nano-clusters shows strong fluctuations, depending on 
cluster size, shape, and position of the atoms in the 
cluster. Nevertheless, trends can be observed, and the 
coupling can be analyzed and understood with respect 
to simple physical terms, namely ferromagnetic double- 
exchange, antiferromagnetic kinetic exchange, hybridiza- 
tion of the d-states, and orientation of the <i-orbitals. As 
these mechanisms act in parallel, they interfere and com- 
pete, giving at the end a complex and seemingly irregular 
picture. It becomes then difficult to isolate the mecha- 
nisms in experiment for a direct testing of the theory. 

We found that the physical picture and the trends be- 
come much more transparent in simple cases, e.g., when 
studying dimcrs or trimers of various configurations. We 
believe therefore that best insight can be gained by study- 
ing highly-symmetric, low-dimensional structures, where 



only one or two parameters define the structure. Such 
could be, for instance, linear vs. corner-shaped trimers, 
or linear vs. zig-zag chains with varying length. As there 
is considerable progress in experimental techniques for 
preparation and probing of such structures, a detailed 
and direct comparison with theoretical results should be 
possible in the near future. 
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